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Abstract. Spectroscopic observations of the chemical abundances in metal-poor H II regions 
provide an independent method for estimating the primordial helium abundance. H II regions 
are described by several physical parameters such as electron density, electron temperature, 
and reddening, in addition to y, the ratio of helium to hydrogen. It had been customary to 
estimate or determine self-consistently these parameters to calculate y. Frequentist analyses 
of the parameter space have been shown to be successful in these parameter determinations, 
and Markov Chain Monte Carlo (MCMC) techniques have proven to be very efficient in 
sampling this parameter space. Nevertheless, accurate determination of the primordial he- 
lium abundance from observations of H II regions is constrained by both systematic and 
statistical uncertainties. In an attempt to better reduce the latter, and continue to better 
characterize the former, we apply MCMC methods to the large dataset recently compiled by 
Izotov, Thuan, &: Stasihska (2007). To improve the reliability of the determination, a high 
quality dataset is needed. In pursuit of this, a variety of cuts are explored. The efficacy of 
the He I A4026 emission line as a constraint on the solutions is first examined, revealing the 
introduction of systematic bias through its absence. As a clear measure of the quality of the 
physical solution, a analysis proves instrumental in the selection of data compatible with 
the theoretical model. In addition, the method also allows us to exclude systems for which 
parameter estimations are statistical outliers. As a result, the final selected dataset gains 
in reliability and exhibits improved consistency. Regression to zero metallicity yields Yp = 
0.2534 lb 0.0083, in broad agreement with the WMAP result. The inclusion of more observa- 
tions shows promise for further reducing the uncertainty, but more high quality spectra are 
required. 
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1 Introduction 

Next to the cosmic microwave background radiation, standard big bang nucleosynthesis 
(SBBN) is the most robust probe of the early universe available [1-3]. Furthermore, us- 
ing the precise baryon density as determined by WMAP [4, 5], SBBN has effectively become 
a parameter free theory [6]. As such, one can use SBBN to make relatively precise predictions 
of the initial light element abundances of D, '^He, ^He, and ^Li [7-16]. Therefore, an observa- 
tional determination of these abundances becomes a test of the concordance between SBBN 
theory and the analyses of microwave background anisotropics. To test these predictions, 
the observed abundances must be determined with relatively high precision. Unfortunately, 
there is a logarithmic relationship between the baryon to photon ratio, rj, and the primordial 
helium abundance, Yp. Thus, any meaningful test of the theory requires a determination 
of Yp to an accuracy of < 1%. The 7-year WMAP value for 77 is (6.19 ± 0.15) x 10~^° [5]. 
For comparison, the SBBN calculation of Cyburt et al. [15], assuming the WMAP r] and a 
neutron mean life of 885.7 it 0.8 s [17], yields Yp = 0.2487 ± 0.0002, a relative uncertainty of 
only 0.08%. 

Here, we discuss the determination of Yp using observations of low metallicity H II re- 
gions in dwarf galaxies. By fitting the helium abundance versus metallicity, one can extrapo- 
late back to very low metallicity, corresponding to the primordial helium abundance [18]. The 
oxygen to hydrogen ratio, 0/H, commonly serves as a proxy for metallicity. The difficulties 
in calculating an accurate and precise measure of the primordial helium abundance are well 
established [19-21]. Previously, we introduced a new method based on Markov Chain Monte 
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Carlo (MCMC) techniques [22, A0S2], and here, we systematically apply this technique to 
the data compiled in [21, ITS07]. 



Observations of the helium to hydrogen emission line ratios from extragalactic H II 
regions provide a measure of the helium to hydrogen ratio, y"*" = -^^j^jjy- The uncertainty 
in y^ follows from the statistical measurement errors of the helium and hydrogen emission 
line fluxes, in addition to a myriad of systematic effects. Interstellar reddening, underlying 
stellar absorption, radiative transfer, and collisional corrections alter the observed flux, com- 
plicating the determination of y"*", and amplifying the uncertainty. First, dust along the line 
of sight scatters the emitted photons (interstellar reddening). Second, the stellar continuum 
juxtaposes absorption features under the emission lines (helium and hydrogen underlying 
absorption). Third, before leaving the H II region the emitted photons can be absorbed and 
re-emitted (radiative transfer). Finally, in addition to the dominant recombination emission, 
collisional excitation also contributes to the emission (collisional corrections for helium and 
hydrogen). Because none of these processes is directly measured, they cannot be removed 
independent of the observed emission lines and theoretical models. As a result, the uncer- 
tainty on y+ must reflect the presence of these systematic effects. Therefore, high fidelity 
spectra are required to accurately determine y"*", while simultaneously estimating the model 
parameters needed to correct for the listed systematic effects. This desired confidence weighs 
against the need for larger sample sizes to decrease the uncertainty on Yp (and dY/dZ). 

Motivated by the work of Izotov, Thuan, & Lipovetsky [23] and Peimbert, Peimbert, & 
Ruiz [24], the importance of Monte Carlo techniques was demonstrated in a "self-consistent" 
analysis method for determining the nebular helium abundance based upon six helium and 
four hydrogen lines [19, 20]. In preceding work. Aver et al. [25, AOS] updated and extended 
the physical model and integrated the helium and hydrogen calculations with the goals of 
increasing accuracy and removing assumptions. Improving the statistical technique of AOS, 
Aver et al. [22, A0S2] used a Markov Chain Monte Carlo (MCMC) analysis to explore a 
global likelihood function for all model parameters including the helium abundance. 

The small sample size employed in AOS & AOS 2 implied a clear avenue for improvement. 
The small sample had limited capacity to constrain the slope of He vs. 0/H and contributed 
to the large uncertainty on the determination of Yp. Since the publication of Olive & Skillman 
[20], Izotov et al. have continued making new observations of H II regions, significantly 
increasing their dataset [21, 26]. This work makes use of that expanded dataset, aided by the 
MCMC analysis, which is used in the process of making quality cuts in the dataset. Included 
in this process is the use of ^ ^ measure of goodness-of-fit between the model and data. 
A new, more rigorously selected dataset holds the promise of an improved determination of 
the primordial helium abundance. 

Section 2 briefly summarizes the model of AOS & A0S2. A broad overview of the 
ITS07 dataset is given in §3, and in §4, the subsets of the observations with and without 
the detection of He I A4026 are compared. Subsequently, further cuts on the dataset are 
investigated. First, providing an indication of the quality of the flt, the of the best-flt 
solution is scrutinized in §5. Second, in §6, ill-behaved points with unreliable determinations 
are exhibited and discussed. Third, §7 examines the distribution of the dataset in terms of 
model parameters used to correct for systematic effects. After making cuts on the dataset, 
Yp is determined in §8. Finally, §9 offers a discussion of the results, their exploration, and 
of further improvements in the determination of the primordial helium abundance. 
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2 Model overview 



This work uses the model introduced in AOS and the MCMC statistical analysis introduced 
in A0S2 and applies them to a larger dataset. The basic definitions are summarized below. 
Please see AOS and A0S2 for full details and discussion. CosmoMC^ is used to perform the 
Markov Chain Monte Carlo analysis: efficiently exploring the parameter space and calculating 
the x^ 

( ^2 

2 ^F(Hf3) F{Hli) 



A 

where the emission line fluxes, F{X), are measured or calculated for six helium lines (A3889, 
4026, 4471, 5876, 6678, and 7065) and three hydrogen lines {Ha, H^, HS) each relative to 
HI3. The ™ ^q. 2.1 runs over all He and H lines and (t(X) is the measured uncertainty 
in the flux ratio at each wavelength. The best-fit solution (minimum x^) is then found and 
frequentist confidence levels are determined from Ax^. 

The calculated He fiux at each wavelength A relative to the flux in HP is given by 



F{H/3) " E{H(3) '^(A)+ai.e(A) '1 + %{H(3) 
The ratio of H fluxes is deflned analogously. 



W{Hl3)+aH{Hf}) ^ ^ 

^' lO-fW^^(^P). (2.3) 



FiX) _ EiX) -'"^n-r 1 + §(A) ,.-/(.)c(..) 



F{HP) E{Hf3) ^(^H«^(A) l + g(F/3) 

The predicted model fluxes shown above are calculated from an input value of y"*" and an 
emissivity ratio of the helium or hydrogen line to H/3, -gp^, with corrections made for red- 
dening (C(H/3)), underlying absorption (an & a^e)) collisional enhancement, and radiative 
transfer. The optical depth function, /,-, and collisional to recombination emission ratio, ^, 
are both temperature (T) and density (ug) dependent (the emissivities are also temperature 
dependent). The parameters a^e and an correspond to A4471 and H/3 respectively. The 
wavelength dependence of the underlying absorption is discussed in detail in AOS. Addition- 
ally, the hydrogen collisional emission depends on the neutral to ionized hydrogen ratio 
Therefore, there are a total of eight model parameters (y"*", Ug, ane, t, T, C(H/?), an, 0- 
extensive description and analysis of the physical model is provided in AOS. The statistical 
method of sampling the multi-dimensional parameter space is described in A0S2. 

The model fluxes also rely on the measured equivalent widths (W{X)). However, the 
flux of the continuum at each wavelength, /i(A), which relates the line flux to the equivalent 
width, is constrained such that changes in the equivalent width are proportional to changes 
in the flux (see A0S2): 

h{X) _ F(A) WiHl3)r 



h{H(3) F{HI3)^^^^ W{X)r 
F{X) W{X) h{X) 



F{H(3) W{HP) h{H(3) 



(2.4) 
(2.5) 
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Finally, the temperature derived from the [O III] AA4363, 4959, and 5007 emission lines, 
T(0 III), is incorporated as a prior with the resulting term added into the y^: 

x'T = iT-T{OIII)f/a^), (2.6) 

where a = 0.2 Tqiii provides a very weak constraint, but is useful in eliminating unphysical 
solutions (see A0S2). 

3 The ITS07 sample 

This work analyzes the HeBCD sample of Izotov, Thuan, & Stasihska [21, ITS07]. The 
HeBCD sample consists of 93 different observations of 86 H II regions in 77 galaxies. As 
described in ITS07, the spectra were reduced and extracted using IRAF. The line flux errors 
are determined from statistical counting uncertainty combined with an assumed standard 
star absolute flux calibration error of 1%. In addition to the helium and hydrogen line flux 
ratios and equivalent widths, this work makes use of T(0 III) and O/H, as reported in 
ITS07. Corrected line fluxes and values for He I A4026 were obtained from Y. Izotov (private 
communication) . 



1 ' 1 ' 1 ' 1 ' r 
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Figure 1. Histogram of the HeBCD sample of ITS07 in terms of the oxygen to hydrogen ratio, 
O/H. 

To begin to characterize the sample, the distribution in O/H is shown in figure 1. 
As lower metallicities are more desirable for determining Yp, the ITS07 sample targeted low 
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Figure 2. Histogram of the HeBCD sample of ITS07 in terms of the equivalent width of H/3, W(H/3) . 

metallicities. Slightly more than half of the sample have 0/H less than the 0/H = 9.2 x 10~^ 
limit imposed in Olive &: Skillman [20]. For reference, the solar value is O / H q= (4.90 ± 
0.02) X 10~^ [27]. This work does not impose a hard cut-off on 0/H. Instead, first the entire 
sample is analyzed, then screened for quality, and the distribution of metallicities is examined 
only as a final step. Ideally, the density of points along the metallicity baseline should be high 
enough to capture the uncertainty in the solutions; as a result, any outliers will be excluded. 

Figure 2 shows the sample's distribution in W(H/3). At lower values of W(H/3) the 
corrections for underlying absorption have a larger effect; as a result, Olive & Skillman 
[20] only used objects with W{H(5) > 200 A. In this work, as with 0/H, all objects are 
analyzed and then investigated for reliability. In particular, the behavior of objects with large 
underlying absorption corrections is examined. Furthermore, given the decreasing signal 
to noise ratio of low equivalent width objects, in addition to their larger susceptibility to 
systematic uncertainty introduced by the corrections for underlying absorption, one could 
imagine identifying a threshold in W(H/3) below which reliable solutions were not possible. 
However, we found no such threshold. Some objects with relatively low equivalent widths 
~ 100) proved to be reliable and showed similar solution characteristics to higher 
W{HI3) spectra. 

4 The importance of He I A4026 

For 23 of the 93 observations in the HeBCD sample, the weak line He I A4026 is not detected. 
These objects can still be analyzed as described in §2, with the helium composed of the sum 
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of the remaining five He lines. To help evaluate the reliability of these 23 objects, the impact 
of removing He I A4026 is investigated in five H H regions. These galaxies were chosen to 
sample a range of equivalent widths: SBS 0917+527, HS 0924+3821, SBS 1152+579, Mrk 209, 
and SBS 0940+544 2, with W(H/3) of 85.9, 114.2, 189.3, 224.1, and 241.1 A, respectively. An 
eight parameter fit to yield the minimum was performed and uncertainties were calculated 
from Ax^. The results of this analysis are given in table 1. Note that without He I A4026, 
the system is, in principle, over-constrained; though due to degeneracies, a "perfect = 0" 
is not found. The absolute and fractional differences between the results, when solved with 
and without He I A4026, are also given in table 1. 

From table 1, we see that the most significant effect of removing He I A4026 is clearly, 
and not surprisingly, on the value of ^uei which increases by 0.23 A, on average, which 
is larger than the typical uncertainty in the solutions including He I A4026. The average 
fractional increase is 69%. The uncertainty on aj^e also increases notably, with the average 
uncertainty tripling. From table 1, it is clear that none of the other determinations of physical 
conditions show a similar significant bias with the absence of a He I A4026 measurement^. 

Since the correction for underlying helium absorption translates linearly into the helium 
abundance, as can be seen in eq. 2.2, increases in the underlying helium absorption directly 
lead to an increased y"*". Table 1 shows an increase in y"*" for all five objects, with an average 
fractional increase of 4%. Clearly this is an unacceptable bias when the ultimate goal is an 
uncertainty on the order of ~ 1%. 

This systematic bias is further demonstrated in figure 3, which shows the helium abun- 
dances for all 93 objects. The set of 23 observations without He I A4026 exhibits bias toward 
higher values of y"*", with a subset of these objects noticeably elevated from the majority of 
the points with He I A4026. A calculation of the mean for each population underscores the 
shift; the mean abundance for the 70 objects with He I A4026 is < >= 0.0867 ± 0.0056, 
while it is < y+ >= 0.0900+0.0352 for the 23 without He I A4026. The uncertainties reported 
with each < y+ > are the dispersions of the sample. The seven-fold increase in dispersion 
for the sample without He I A4026 demonstrates its unreliability. 

He I A4026, being the weakest of the six helium lines used in our analysis, is the most 
sensitive to underlying helium absorption. He I A4026 also shows relatively weak sensitivity 
to temperature, density (collisional enhancement), and optical depth, and Olive &: Skillman 
[19] demonstrated its usefulness for the purpose of constraining underlying helium absorption. 
However, this also requires that the input spectrum is a very high quality one. The addition 
of even weaker He I lines such as A4922, A7281, and A4387 [e.g., 24, 28], which are all singlet 
lines with low susceptibility to optical depth effects, could, in principle, provide even stronger 
constraints on the effects of underlying absorption. For this to be true, the spectra would 
need to be of much higher quality than is typically obtained for this type of work. 



^While aH also shows a systematic increase when A4026 is removed, the shift is small compared with the 
original uncertainties. 
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Table 1. Comparison of the solutions for five sample galaxies with and without the inclusion of the He I A4026 observation. Units used in the table 
are: Ug (cm"'^), a.He and a,H (A), Te (K). Other quantities are dimensionless. 
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Figure 3. Plot of the helium abundance versus 0/H for the entire sample. The black circles denote 
objects for which He I A4026 is reported and was used to calculate and determine the best-fit 
solution, while the red squares mark objects for which He I A4026 is not reported and is, therefore, 
excluded from their analysis. An upward bias in the value of y"*" is apparent in the distribution of 
points lacking He I A4026. 



The analyses above indicate the importance of including the faint He I A4026 emission 
line in analyses like ours where the underlying absorption is solved for in a minimization. 
Solutions lacking it are prone to drifting to larger values of a//e and, consequently, y"*". As 
a result of the testing above, and the corresponding behavior manifested in the population 
of objects without He I A4026, these 23 objects without He I A4026 detections are excluded 
from the subsequent analysis and results. A primary goal of this work is to reduce systematic 
bias as much as possible, and the strongest constraint upon a/^e is provided by He I A4026. 
As is evident in table 1, the benefit of an accurate He I A4026 measurement to determining 
the helium abundance should not be underestimated, and thus, requiring it improves the 
reliability of the sample. Note that the lack of A4026 in other analyses does not automatically 
result in artificially higher values of y"*"; however, not properly accounting for underlying 
helium absorption will certainly affect the reliability of those results. 

5 Using chi-squared as an analysis tool 

After dropping the 23 objects without A4026 measurements, 70 objects remain. Figure 4 
shows the distribution of values for the solutions for these remaining 70 objects. Clearly, 
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not all solutions provide good fits to the data (as evidenced by relatively high values of x^)- 
Fits with large may indicate a measurement discrepancy in the line fluxes, an underesti- 
mation of the uncertainties, or the possibility that the model used to derive abundances is 
inappropriate for the object under study. In either case, points with large are not reliable, 
and thus, we exclude them from further analysis in determining Yp. There are nine observed 
line ratios used to calculate x^ and eight model parameters fit to the data; thus, there is 
only one degree of freedom, modulo correlations^. Here we choose to cut the sample at a 
standard value of x^ < 4, corresponding to a 95.45% confidence level^. This cut removes 45 
observations, leaving 25 (figure 4). The disappointing result that nearly two-thirds of the 
sample have solutions with such low likelihoods (x^ > 4) is troubling and warrants further 
investigation. At this point, we cannot be certain whether this result is due to deficiencies 
in the model, the observations and their errors, or both. 

Since the number of objects with high values of x^ is so large, in the following sections, 
we will carry out parallel analyses in which we study the reduced sample of 25 objects and 
compare to the sample of 70. What we find is that a majority of the solutions with high x^ 
values also tend to have questionable values for different physical parameters. In the end, we 
conclude that cutting the sample on x^ is fully justified on both theoretical and empirical 
grounds. 

6 Investigating problematic points 

Next, we further investigate the quality of the solutions for the 25 objects meeting the 95% 
confidence level. Clearly it makes sense to exclude non-physical or ambiguous solutions even 
if those solutions have acceptable values of x^- As was demonstrated in A0S2, one of the 
benefits of an MCMC analysis is that likelihood plots illuminate deficiencies in an object's 
solution, such as very weakly constrained parameters and multiple minima. For example, the 
neutral to ionized hydrogen fraction, ^ is of particular concern. As was discovered in AOS, 
and further investigated in A0S2, ^ is very weakly constrained at relatively low temperatures 
(higher metallicity) , and can admit completely unphysical solutions for the relative amount 
of neutral hydrogen. 

Figure 5 illustrates the difference between the well-constrained, physically realistic so- 
lution of Mrk 209, and the unbounded solution of Mrk 35. For Mrk 35, the global minimum 
has not yet been reached at a value of ^ = 100, which corresponds to 99% of the hydrogen 
being neutral (^ = lO"'^ is a characteristic value for H II regions). The best fit solutions 
with high neutral fraction are not consistent with classifying this object as an H II region (or 
photo- ionized) ! Observations that do not yield a physically meaningful solution imply either 
model deficiencies or errors in the observed spectrum. Thus, we implement a very conser- 
vative criterion that any object whose solution exceeds ^ = 0.333 (25% neutral hydrogen) is 
excluded from further analysis. There are two objects from the reduced sample of 25 with 
large neutral fractions that we exclude with this cut, and thus, we go forward with a reduced 
sample of 23. 

Figure 6 shows the distribution of values of ^ in the original sample of 70, prior to the 
X^ cut. Note that there are 19 solutions with ^ > 0.333 in this original sample, 17 of which 
have x^ > 4. Overall, unrealistically large values of the neutral fraction are indicative of low 

^Correlations between the parameters will increase the effective degrees of freedom. 

*We note that of the 23 objects lacking A4026, nine had > 4 (though, as noted, in the absence of 
degeneracies among the parameters, we should expect a of by fitting 8 parameters to 8 observables.) 
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Figure 4. Histogram of best-fit for the 70 observations with He I A4026 reported. 25 observations 
have < 4 (95.45% confidence level). 



probability solutions. We also identify objects with ^ > 0.01, and whose lower error bound 
does not encompass = 0.001. These solutions are flagged for special special consideration 
and will be discussed later. There are three such objects in the reduced sample of 23 (seven 
before the cut). 

Additionally, the likelihood plots of each object are inspected for ambiguous solutions. 
Only two objects yielding uncertain solutions are found, both of which have very large x^- 
SBS 1420+544 exhibits a pronounced double minimum, the lower of which corresponds to a 
very low temperature (~10,000 K). Such behavior calls into question the reliability of this 
object, and is particularly striking as it is the only object in the dataset to exhibit this phe- 
nomenon. Examples of solutions with double minima were discussed in A0S2. J0519+0007 
is the only object to exhibit a likelihood for y"*" which is not strongly parabolic about its 
minimum. Shown in figure 7, the abundance is nearly unbounded for values greater than the 
minimum, and in fact, the values of y"*" attained are larger than for any of the other objects. 
Mrk 209 is again provided for comparison to a normal, well-behaved abundance. The values 
of for SBS 1420+544 and J0519+0007 - 15.7 (for the more physical minimum) and 18.3, 
respectively ~ provide a strong sign of their unreliability. Similarly, it is interesting that both 
objects have large optical depth values (4.3 & 5.6). Furthermore, J0519+0007 has a temper- 
ature much lower than its 0[III] temperature (12,800 K vs. 20,700 K) and an unphysically 
large neutral hydrogen fraction (,^=0.67). 
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Figure 5. versus the neutral to ionized liydrogen fraction, ^, for Mrk 35 and Mrk 209. On tire 
left, Mrk 35 exhibits an unphysical ^ in contrast to the well determined value for Mrk 209, shown on 
the right, with the 68% confidence level marked. Mrk 35's ^ likelihood plot corresponds to more than 
99% of the hydrogen present being in the neutral state, an unphysical model for an H II region. 
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Figure 6. Histogram of the neutral hydrogen to total hydrogen fraction, ^/(l + ^), for the 70 
observations with He I A4026. The 25 with < 4 are over-plotted using darker, cross-hatched bars. 



After removing the objects with problematic solutions we have two samples. There are 
50 objects with well defined solutions, and 23 of these have satisfactorily low values of x^- 
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Figure 7. versus helium abundance (y+) for J0519+0007 and Mrk 209. On the left, J0519+0007 
exhibits a very weakly bounded y^ in contrast to the well determined value for Mrk 209, shown on the 
right, with the 68% confidence level marked. J0519+0007's y+ likelihood plot is symptomatic of an 
object whose determination fails on multiple different physical grounds and is completely unreliable. 



7 Looking for model outliers 

With a dataset comprised of 23 observations, for which the model is a good fit and the 
parameters are clearly determined, Yp can be extracted. However, the models used for 
correcting for observed flux for optical depth and underlying absorption carry significant 
systematic uncertainties. The effect of these systematic uncertainties can be minimized by 
limiting the size of the corrections. Additionally, the solution temperature and the 0[III] 
temperature should be in relatively good agreement for the solution to make sense. Here we 
search for anomalous values of all three physical parameters. We also discuss the choice of a 
metallicity baseline. 

7.1 Anomalously large values of optical depth 

In A0S2, it was noted that high values of optical depth increase susceptibility to model 
deficiencies. Namely, the radiative transfer calculations of Benjamin, Skillman, & Smits [29] 
do not take into account any expansion, non-uniformity, or turbulence in the H 11 region. A 
histogram showing the distribution of the optical depth in objects before and after the < 4 
cut is shown in figure 8. There are originally 3 potential outliers with r > 4. Reassuringly 
though, the cut removes two of these. The remaining object is flagged to ascertain the 
impact of questionable points on the regression. 

7.2 Anomalously large values of underlying absorption 

Next, a similar investigation is performed for underlying absorption. Though AOS added 
wavelength dependence to the modeling of underlying absorption, it is still a difficult effect 
to accurately quantify. It was for this reason that Olive &: Skillman [20] only analyzed galaxies 
with W{HI3) > 200 A, which minimizes the effect of underlying absorption. The histograms 
for an and ajje are shown in figures 9 and 10, respectively. As one can see, there are four 
potential outliers with an > 6 A and two with ajje > 1 A. One of the objects is anomalously 
high in both, and one is removed by the selection. Thus, four objects are flagged for 
further analysis with acceptable values of but anomalous solution parameters. 
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7.3 Anomalous values of temperature 



The comparison of the solution temperature and the 0[III] temperature provides a final 
physical parameter to investigate here. In this case, the issue is not one of outliers and their 
model realizability, but the physical constraint on these two temperatures: they should be 
close to one another, with T(0 III) serving as a loose upper bound for T. Note that T is 
very weakly constrained by T(0 III), but primarily determined by the helium lines (i.e., 
T~T(IIe II); see A0S2 for further discussion). 

Figure 11 shows the difference, AT = T{0 III) — T, before and after the cut. Since, 
in the presence of temperature fluctuations, the 0[III] temperature is biased toward higher 
values, T(He II) is estimated to be less than T(0 III) by 3-11% (see A0S2). Additionally, 
the uncertainties on T are very large, averaging ~2000 K over the set of 50 observations. As 
a result, there will be a spread in AT, including negative values. However, the most extreme 
values for this difference, AT < —3000 K and AT > 5000 K are very unlikely results (as is 
the number of results with AT < K). Reinforcing the utility of the selection, the < 4 
subset does not admit the extreme values, and the distribution now reflects the physical 
expectation. 
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7.4 Choosing a baseline in metallicity 



The last criterion for dataset selection is the metallicity baseline. Extending the metallicity 
baseline to higher values increases the systematic uncertainty in the value of the primordial 
helium abundance which is due to the assumption of a linear relationship of He/H with 
chemical evolution. For this reason, a shorter metallicity baseline is preferred, but there is 
a trade-off between a metallicity cut and sample size. Clearly it is not desirable to truncate 
the dataset significantly and thus lose precision on the determination of Yp. Figure 12 shows 
that both samples exhibit a reasonable coverage of points up through O/H = 15.1 x 10~^ 
and then a possible outlier at 18.1 x 10~^. This last point could have undue leverage on the 
regression if included and is therefore excluded. The exclusion of the largest metallicity point 
yields a final dataset of 22 objects with well-defined solutions that are good fits to the data. 

A summary of the cuts and flagged objects on the sample is presented in table 2. To 
emphasize the utility of the goodness-of-fit test, the sequence is provided with < 4 
applied first, as was presented in the above description, and after the other criteria. The 
resulting final dataset is the same, but the effectiveness of for identifying problematic 
solutions is apparent. 
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Figure 10. Histogram of tlie underlying hydrogen absorption, a/fe, for tlie 50 well-defined solutions. 
The 23 with < 4 are over-plotted using darker, cross-hatched bars. 

8 Results from the Final Dataset 

The 22 objects for which the model is a good fit, which return physically meaningful param- 
eter solutions, and which provide a robust metallicity baseline for regression comprise the 
Final Dataset, which we use to determine Yp. Each of the objects, with the results of their 
best-fit solutions and uncertainties, are presented in table 3. Figure 13 presents the derived 
y+ values as a function of 0/H. The seven objects flagged for large outlier values of r, ajy, 
aj^e, and ^ are highlighted with different symbols in figure 13. These flagged data points 
show a possible systematic shift to larger values of y'^; the average value of y"*" for the flagged 
objects is 9% higher than that of the unflagged objects in the Final Dataset. 
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Figure 11. Histogram of tlic difference in the 0[III] temperature and solution temperature, 
T4(0 ///) - T4, wliere T4 is in units of lO^'K, for the 50 weU-defined solutions. The 23 with < 4 
are over-plotted using darker, cross-hatched bars. It is expected that T(0 III) should be greater than 
T, and that T should closely track T(0 III). That some values of T wih be larger than T(0 III) is 
a result of the combination of the close relationship of their actual physical values and the relatively 
large uncertainty on T. However, the most extreme values are not easily justified. As result, it is not 
surprising that the cut removes most of them. 
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Figure 12. Histogram of the oxygen to hydrogen ratio, 0/H, for the 50 well-defined solutions. The 
23 with < 4 are over-plotted using darker, cross-hatched bars. 



-17- 



Criterion 


% < 4 iirst 


It X < 4 reserved 






until last 




yd 


93 


Cuts 






He I A4026 Not Detected 


23/93 


23/93 


2 ^ A 

X > 4 


45/70 




^ > u.ooo 


Z/ ZD 


iy/ /u 


Degenerate Solution 


0/23 


1/51 


Subtotal: Well-Defined Solutions 


23 


50 


Flagged 






C > 0.01 


3 


7 


T > 4 


1 


3 


an > 6 A 


3 


4 


a/Ze > 1 A 


2 


2 


T{0 III) - T < -3000 K 





1 


T{0 III) - T > 5000 K 





1 


Subtotal: Unique Flagged 


8 


17 


Cuts 






0/H > 15.2 X 10-5 


1/23 


1/50 


> 4 




27/49 


Final Dataset 


22 


22 


Flagged 


8 


8 


Qualifying 


14 


14 



Table 2. Breakdown of the cuts and flags on the sample. To better highlight the behavior of the 
physical parameters and emphasize the role of x^, each criterion is tabulated for the ease where the 
cut on is made first and the case where it is made afterward. 
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Object He+/H+ r Te C(H^) &h C x 10 X 



Final Dataset Not Flagged 



I Zw 18SE 1 


0, 


.0811 


+U.UU4ti 
-0.0068 


4 +^48 


0, 


.29 


+U.23 
-0.24 


0, 


.43 


+U.VO 
-0.43 


18,190 


+2500 
-2960 


0, 


.01 


+U.U2 
-0.01 


3, 


.80 


+U.B2 
-0.69 





+26 
-0 


0.5 


SBS 0940+544 2 


0, 


.0870 


+0.0057 
-0.0050 


8+r 


0, 


.51 


+0.18 
-0.18 


0, 


.08 


+0.64 
-0.08 


18,380 


+ 1330 
-2500 


0, 


.07 


+0.03 
-0.05 


1. 


.81 


+ 1.34 
-1.09 


8 


+74 
-8 


0.7 


Tol 65 


0, 


.0825 


+0.0052 
-0.0038 


211 ^ 


0, 


.67 


+0.11 
-0.08 


0, 


.00 


+0.47 
-0.00 


17,330 


+ 1620 
-2910 


0, 


.13 


+0.02 
-0.05 


3. 


.54 


+0.82 
-0.60 


1 : 


+ 137 
-1 


2.4 


SBS 1415+437 (No. 2) 


0, 


.0861 


+0.0029 
-0.0063 


r +343 
-5 


0, 


.54 


+0.07 
-0.12 


1, 


.16 


+0.75 
-1.16 


15,440 


+ 1910 
-2320 


0, 


.00 


+0.02 
-0.00 


3, 


.48 


+0.72 
-0.72 





+57 
-0 


1.9 


HS 1442+4250 


0, 


.0830 


+0.0064 
-0.0068 


1 +257 
^ -1 


0, 


.46 


+0.40 
-0.38 


0, 


.29 


+0.85 
-0.29 


14,880 


+2180 
-2650 


0, 


.03 


+0.03 
-0.03 


0, 


.00 


+ 1.15 
-0.00 


127 


+ 1550 
-127 


2.2 


Mrk 209 


0, 


.0848 


+0.0022 
-0.0052 


I +293 


0, 


.37 


+0.11 
-0.16 


0, 


.49 


+0.41 
-0.49 


16,030 


+ 1380 
-2580 


0, 


.01 


+0.02 
-0.01 


2. 


.16 


+0.95 
-0.80 





+59 
-0 


1.4 


Mrk 71 (No. 1) 


0, 


.0871 


+0.0058 
-0.0057 


93+i« 


0. 


.60 


+0.28 
-0.25 


1, 


.75 


+0.67 
-0.83 


14,170 


+ 1860 
-2090 


0, 


.09 


+0.04 
-0.05 


1. 


.13 


+ 1.71 
-1.13 


132 


+ 1270 
-132 


1.2 


SBS 0917+527 


0, 


.0866 


+0.0069 
-0.0067 


1 +196 
^ -1 


0. 


.16 


+0.13 
-0.12 


0, 


.00 


+0.59 
-0.00 


12,770 


+ 1450 
-1830 


0, 


.05 


+0.05 
-0.05 


0. 


.45 


+0.65 
-0.45 


711 


+5510 
-711 


1.2 


SBS 1152+579 


0, 


.0891 


+0.0066 
-0.0056 


1 +130 
^ -1 


0. 


.42 


+0.15 
-0.14 


2, 


.48 


+0.75 
-0.61 


15,170 


+ 1840 
-2040 


0, 


.23 


+0.04 
-0.05 


0, 


.44 


+0.81 

-0.44 


77 


+498 
-77 


3.8 


SBS 1054+365 


0, 


.0921 


+0.0038 
-0.0085 


1 +428 
^ -1 


0, 


.49 


+0.18 
-0.24 


0, 


.51 


+0.57 
-0.51 


12,040 


+ 1720 
-1620 


0, 


.00 


+0.07 
-0.00 


2. 


.87 


+0.49 
-1.04 


2040 


+5160 
-2040 


1.0 


SBS 0926+606 


0, 


.0966 


+0.0038 
-0.0099 


+296 


0, 


.70 


+0.13 
-0.20 


0, 


.12 


+0.61 
-0.12 


12,790 


+ 1520 
-1750 


0, 


.00 


+0.07 
-0.00 


0, 


.75 


+0.33 
-0.75 


1460 


+4470 
-1460 


0.7 


SBS 1135+581 


0, 


.0854 


+0.0038 
-0.0034 


468"+^^^ 


0. 


.44 


+0.09 
-0.08 


0, 


.00 


+0.59 
-0.00 


11,270 


+ 1880 
-1660 


0, 


.10 


+0.02 
-0.03 


3, 


.42 


+0.47 

-0.54 


^ 


-3160 

-0 


2.7 


HS 0924+3821 


0, 


.0851 


+0.0052 
-0.0049 


QO +836 
»^ -92 


0. 


.35 


+0.13 
-0.12 


0, 


.21 


+0.82 
-0.21 


11,600 


+ 1480 
-2190 


0, 


.16 


+0.02 
-0.05 


2. 


.29 


+0.71 
-0.52 


ol 


10100 



0.9 


UM 439 


0, 


.0906 


+0.0125 
-0.0078 


185 


0, 


.01 


+0.24 
-0.01 


3 


.29 


+0.95 
-1.01 


14,580 


+2040 
-3290 


0, 


.20 


+0.05 
-0.07 


0. 


.57 


+ 1.16 
-0.57 


218 


+2980 
-218 


2.9 



Final Dataset with Flags 



SBS 0335-052E1 


0, 


.0863 


+0.UU8U 
-0.0081 


4+182 





.24 


+U.18 

-0.20 


4. 


.96 


+0.89 
-1.01 
+0.50 

-0.55 


18,120 


+2330 
-2680 


0, 


.09 


+U.Ub 
-0.06 


1.51 


+ 1.21 

-1.25 


17 


+ 126 
-17 


0.9 


HS 0122+0743 


0, 


.0996 


+0.0072 
-0.0092 


1 +97 
^ -1 


1 


.40 


+0.21 
-0.23 


0. 


,55 


18,070 


+2270 
-1790 


0, 


.06 


+0.06 
-0.06 


4.34 


+1.47 

-1.57 


38 


+ 109 
-38 


1.1 


UM 461 


0, 


.0932 


+0.0043 
-0.0120 


1 +295 
^ -1 





.90 


+0.27 

-0.35 


2. 


,07 


+0.67 
-0.97 
+0.52 
-0.99 


18,340 


+2110 
-2850 


0, 


.00 


+0.07 
-0.00 


7.72 


+0.87 
-1.94 


28 


+98 
-28 


1.5 


HS 0811+4913 


0, 


.1004 


+0.0042 
-0.0117 


+323 
^" -26 


1 


.39 


+0.25 
-0.46 


0. 


,99 


14,850 


+ 1750 
-2070 


0. 


.01 


+0.07 
-0.01 


9.13 


+0.84 
-2.89 


396 


+ 1370 
-396 


1.5 


UM 238 


0, 


.0856 


+0.0092 
-0.0060 


1150 i^™ 


0, 


.92 


+0.62 
-0.46 


0. 


,53 


+0.98 
-0.53 


12911 


+ 1950 
-2770 


0, 


.24 


+0.04 
-0.08 


13.21 


+3.69 
-2.19 


178 


+4280 
-178 


0.9 


UGC 4483 


0, 


.0964 


+0.0053 
-0.0041 


1 +113 
^ -1 


0, 


.47 


+0.09 
-0.08 


0. 


,36 


+0.56 
-0.36 


14,090 


+ 1600 
-1680 


0. 


.09 


+0.03 
-0.04 


0.00 


+0.54 
-0.00 


716 


+2070 
-525 


2.2 


SBS 1331+493 


0, 


.0906 


+0.0051 
-0.0065 


iqo +507 
i'iO -138 


0, 


.23 


+0.10 
-0.13 


0. 


,97 


+0.58 
-0.97 


12,910 


+ 1940 
-1570 


0, 


.00 


+0.04 
-0.00 


0.29 


+0.76 
-0.29 


1440 


+3720 
-1170 


3.8 


SBS 1533+574B 


0, 


.0957 


+0.0203 
-0.0133 


401 +^2* 

^Ul _4oi 





.46 


+0.20 

-0.20 


1. 


,53 


+ 1.44 

-1.32 


12,620 


+2750 
-1460 


0, 


.27 


+0.11 

-0.10 


2.56 


+0.95 

-1.03 


2630 


+3690 
-2450 


2.7 



Table 3. Physical Conditions and He+/H+ Abundance Solutions of Final Dataset 
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Figure 13. Plot of y+ vs 0/H for the 22 objects meeting the prescribed rehabihty standards. 
The upward triangles signify points flagged for large outlier values in optical depth or underlying 
absorption. The downward triangles signify points flagged for large neutral hydrogen fractions. 



The primary goal of this work, the primordial helium abundance (mass fraction), Yp, 
can now be calculated for several subsets of the final dataset. A regression of Y, the helium 
mass fraction, versus 0/H, the oxygen to hydrogen mass fraction, is used to extrapolate to 
the primordial value^. The 0/H values are taken directly from ITS07. 

Because it minimizes confounding systematic effects, our preferred dataset is the 14 
qualifying points. The relevant values for its regression are given in table 4. The regression 
yields, 

= 0.2534 ± 0.0083, (8.1) 

with a slope of 54 it 102 and a of 2.9. The result is shown in figure 14. Note that the 
expected value of foi" this dataset is ~ 12, so the resultant is unexpectedly low. This 
result for Yp agrees well with the WMAP value of Yp = 0.2487 ± 0.0002. A0S2 determined 
Yp = 0.2609±0.0117 (0.2573;o:oog8 with the slope restricted to be positive). Given their large 
uncertainties, these results are in agreement with the newer result. The smaller uncertainty 
on the unconstrained fit is a direct result of the increased sample size. 

As the 0/H domain is limited, an estimate of Yp using the mean value is justified and 
gives, 

Yp = 0.2574 ± 0.0036. (8.2) 

^This work takes Z = 20{O/H) such that Y = 4y(i-20(o/g)) 
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Figure 14. Helium abundance (mass fraction) versus oxygen to hydrogen ratio regression calculating 
the primordial helium abundance. 



This is not significantly different from the result of the regression fit; however, the uncertainty 
is decreased by more than a factor of two. 

Including the flagged objects raises the intercept and reduces error to 0.2611 ± 0.0067 
with a slope of it 86. The reduced uncertainty is a result of the increased number of 
points in the regression, and the possible systematic bias toward larger y"*" within the flagged 
dataset raises the intercept. Olive & Skillman [20] restricted the metallicity baseline to 
O/H = 9.2 X 10^^ Adopting the same metallicity cut with the dataset of this work decreases 
the intercept substantially to 0.2465 it 0.0134 and produces a strongly positive slope (though 
still consistent with 0) of 196 it 230. Using all 93 observations included in their HeBCD 
sample, ITS07 determined Yp = 0.2516 =b 0.0011. Their much smaller uncertainty is achieved 
primarily though the use of the full sample of observations. Table 5 summarizes the calculated 
regression Yp and slope, as well as the mean, <Y>, for several subsets of the Final Dataset 
found in this work. 



9 Summary 

Because of its importance as a probe of big bang cosmology, the predictions of BEN are 
under constant scrutiny from new observations. With the baryon-to-photon ratio fixed by 
the WMAP determination of the baryon density, BEN leads to a distinct set of predictions 
for the light element abundances. The prediction of the primordial D/H abundance is a 
major success when compared to determinations of D/H from quasar absorption system 
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Object 


He+ /H+ 


He++/H+ 


Y 


0/H X 10^ 


I Zw 18SE 1 


0.08112 ± 0.00680 


0.0008 ± 0.0008 


0.2467 ± 0.0155 


1.5 ± 0.1 


SBS 0940+544 2 


0.08701 ± 0.00567 


0.0005 ± 0.0005 


0.2592 ± 0.0125 


3.2 ± 0.1 


Tol 65 


0.08247 ± 0.00520 


0.0011 ± 0.0011 


0.2503 ± 0.0119 


3.5 ± 0.1 


SBS 1415+437 (No. 2) 


0.08610 ± 0.00625 


0.0000 


0.2560 ± 0.0138 


4.2 ± 0.1 


HS 1442+4250 


0.08295 ± 0.00676 


0.0026 ± 0.0026 


0.2547 ± 0.0161 


4.4 ± 0.2 


Mrk 209 


0.08475 ± 0.00518 


0.0010 ± 0.0010 


0.2552 ± 0.0117 


6.1 ± 0.1 


Mrk 71 (No. 1) 


0.08713 ± 0.00580 


0.0008 ± 0.0008 


0.2598 ± 0.0128 


7.2 ± 0.2 


SBS 0917+527 


0.08664 ± 0.00686 


0.0020 ± 0.0020 


0.2613 ± 0.0155 


7.4 ± 0.3 


SBS 1152+579 


0.08906 ± 0.00659 


0.0012 ± 0.0012 


0.2648 ± 0.0144 


7.7 ± 0.2 


SBS 1054+365 


0.09214 ± 0.00850 


0.0000 


0.2688 ± 0.0181 


9.6 ± 0.4 


SBS 0926+606 


0.09660 ± 0.00987 


0.0008 ± 0.0008 


0.2799 ± 0.0205 


10.0 ± 0.4 


SBS 1135+581 


0.08538 ± 0.00382 


0.0008 ± 0.0008 


0.2558 ± 0.0086 


11.7 ± 0.3 


HS 0924+3821 


0.08509 ± 0.00516 


0.0000 


0.2533 ± 0.0115 


12.3 ± 0.5 


UM 439 


0.09064 ± 0.01247 


0.0000 


0.2654 ± 0.0268 


12.3 ± 0.3 



Table 4. Primordial Helium Regression Values 



Dataset N Yp dY/d{0/H) <Y > 

Qualifying 14 0.2534 ± 0.0083 54 ± 102 0.2574 ± 0.0036 

Qualifying + Flagged, excluding ^ > 0.01 19 0.2580 ± 0.0074 16 ± 90 0.2592 ± 0.0034 

Qualifying + ^ > 0.01 Flagged only 17 0.2598 ± 0.0074 + 95 0.2598 ± 0.0032 

Qualifying + All Flagged 22 0.2611 ± 0.0067 + 86 0.2611 + 0.0031 

Qualifying, O/H < 9.2 x 10"^ (A0S/A0S2) 9 0.2465 ± 0.0134 196 + 230 0.2564 + 0.0045 



Table 5. Comparison of Yp for selected datasets 



observations. To further test the theory, accurate '^He abundances are necessary as the 
BBN prediction for ^He is at the 0.1% level. Selecting an improved dataset from existing 
observations was a primary aim of this work. 

Unfortunately, ^He abundance determinations from emission line fluxes carry significant 
systematic uncertainties. In A0S2, we showed that one's ability to accurately estimate the set 
of physical parameters which are used to determine the '^He abundance is greatly improved 
with MCMC methods. An eight-dimensional parameter space (which includes the helium 
abundance) is sampled to calculate a set of line fluxes which can be compared directly with 
observed fluxes. 

Here, we have applied the MCMC method to a large sample of 93 objects, and then 
carefully screened these objects to remove systematic bias, poor quality model fits, and 
spurious physical results. In this process, the importance of He I A4026 was once again 
demonstrated. Data lacking this line was found to be susceptible to a systematic increase 
in underlying Helium absorption and a corresponding bias to larger abundances. 70 of 93 
objects in the initial sample contain the A4026 line, and were the subject of our further 
analysis. 

One of the benefits of the MCMC method is our ability to directly test the goodness- 
of-fit for a particular set of solutions. While we are able to find a best fit solution almost 
every object, there is no guarantee that that solution is actually a good (in a statistical 
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sense) description of the data. As the best fit solution (fitting nine observables with eight 
parameters) is found by minimizing x^, the resulting minimum is a direct measure of the 
goodness-of-fit. Unfortunately, 45 of the 70 objects resulted in a per degree of freedom > 4. 
In our quest for high quality and reliable ^He data, these data were excluded from further 
analysis. On a positive note, the practical impact of the 95% confidence level cut preferentially 
removed outlier parameter solutions. Indeed, only a handful of objects surviving the were 
flagged as outliers. 

The cumulative effect of the careful selection led to a dataset with 14 objects. However, 
the reliability of this Final Dataset has been methodically evaluated and is a distinct im- 
provement over previous efforts. The precision of the resulting primordial helium abundance 
is increased, even if not to the level we hope to someday achieve. The concordance of the 
various regression datasets incorporating fiagged points underscores the robustness of the 
result and the improved quality of the dataset. 

As we have emphasized before, the most promising avenue for significant future im- 
provement in the primordial helium abundance determination lies with higher quality, high 
resolution spectra. Decreased measurement uncertainties will clearly better constrain the 
solution and decrease systematic uncertainty. Beyond that straightforward benefit, there are 
a couple of promising, related paths. First, as discussed in AOS, high resolution spectra 
afford the chance to measure the absorption underlying the Balmer lines directly. Removing 
underlying hydrogen absorption as a solution model parameter would reduce degeneracies 
and produce better constrained parameters. Similarly, high signal to noise spectra open up 
the possibility of adding a weaker Helium or Hydrogen line to our analysis. Again, increas- 
ing the number of lines constraining the parameters would reduce degeneracies and decrease 
systematic uncertainty. That both of these improvements would directly benefit the neutral 
hydrogen fraction, the least well constrained parameter, provides additional impetus. Finally, 
we also plan to quantify the correlations between the parameters for each object. This will 
yield further insight into the goodness-of-fit. 

In summary, we have demonstrated the rigor and transparency of MCMC methods 
in selecting the best available data to ultimately extract the primordial He abundance. It 
supports a stringent screening of candidate spectra, yielding a robust sample. The fruit 
of these labors is an improved determination of the primordial helium abundance, both 
in increased confidence in its accuracy and in a modest increase in its precision. These 
uncertainties are still relatively large, however, and the case for needing higher quality spectra 
is further strengthened. 
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